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We study the driven translocation of polymers under time-dependent driving forces using N- 
particle Langevin dynamics simulations. We consider the force to be either sinusoidally oscillating 
in time or dichotomic noise with exponential correlation time, to mimic both plausible experimental 
setups and naturally occurring biological conditions. In addition, we consider both the case of 
purely repulsive polymer-pore interactions and the case with additional attractive polymer-pore 
interactions, typically occurring inside biological pores. We find that the nature of the interaction 
fundamentally affects the translocation dynamics. For the non-attractive pore, the translocation 
time crosses over to a fast translocation regime as the frequency of the driving force decreases. In 
the attractive pore case, because of a free energy well induced inside the pore, the translocation 
time can be a minimum at the optimal frequency of the force, the so-called resonant activation. In 
the latter case, we examine the effect of various physical parameters on the resonant activation, and 
explain our observations using simple theoretical arguments. 

I. INTRODUCTION 

Translocation of polymers across a nanopore is a ubiquitous process in biology, with examples such as DNA and 
RNA transport through nuclear pore complex, protein transport through membrane channels, and virus injection 
into cells pQ. Kasianowicz et al. [2] demonstrated in vitro that an electric field can transport single-stranded (ss) 
nucleotides through an a-hemolysin membrane channel and it is possible to characterize individual molecules by 
measuring the ionic current blockade when the chain moves through the pore. Later, Li et al. showed [3] that also 
solid-state nanopores can be used for similar experiments with a tunable size of the pore. To further the understanding 
of numerous biological processes and examine the perspective of technological applications such as sequencing and 
gene therapy, there have been extensive experimental [H-TiT] and theoretical studies [P2Tl37| . 

One of the most important quantities of the process is the translocation time and its dependence on the various 
system parameters such as chain length, type of driving force, pore width, etc. Even with the same chain lengths, 
recent experiments [IHZj have shown that different nucleotides exhibit unique patterns in, e.g., the translocation time 
distribution. In particular, Meller et al. [SJ [7] have shown that in the translocation can discriminate between poly- 
deoxyadenylic acid (poly(dA)) and polydeoxycytidylic acid (poly(dC)) with the same chain length. The translocation 
time of poly(dA) is found to be longer with an exponential distribution while that of poly(dC) is shorter with a narrow 
distribution. The origin of the different behavior for each nucleotide was attributed to different interaction between 
the polymer and the pore. Recent simulation studies of Luo et al. [22, 25 quantitatively support this idea. 

Until now, most of the in vitro experimental as well as theoretical studies of polymer translocation have used 
static driving forces. However, it could be important to consider time-dependent forces to understand the process 
in vivo. In a cellular environment the driving forces can be time-dependent due to the nonequilibrium fluctuations 
in the membrane potential, fluid density, and ionic strength, etc. In the case of translocation driven by a molecular 
motor [38], depending on the ATP concentration the driving force can also fluctuate. Motivated by these facts, Park 
and Sung [3j3] studied the translocation of a rigid rod in the presence of a dichotomically fluctuating force. They found 
that the system exhibits resonant activation [40] . where the translocation time attains a minimum at an optimum 
flipping rate of the dichotomic force that is comparable to the translocation rate in the absence of the force. Although 
the study gives valuable insight on the effects of fluctuating forces in polymer translocation, the study is somewhat 
limited, however, as the flexibility of the chain is not considered and, the reflecting boundary condition which forbids 
the chain escape to the cis side is in many cases artificial. Recent molecular dynamics simulations [41] have shown 
that an alternating electric field in a nanopore exhibits a unique hysteresis in the nucleotide's dipole moment and 
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in the chains back-and-forth motion arising from the reorientation of the DNA bases in the nanopore constriction. 
The authors suggest detection of DNA sequences by measuring the potential or the change of the DNA mobility in 
the pore. This study indicates that a time-dependent driving force may be useful for technological application as 
well. In addition, recent Langevin dynamics study shows that the translocation time can be significantly shortened 
by oscillations of the pore width |42| , which accentuates the importance of polymer-pore interactions in the problem. 

In addition to polymer translocation, there are a few simulation studies of different types of polymer transport 
in the presence of time-dependent driving forces. Tessier and Slater [43] considered polymer transport through a 
microchannel in the presence of a periodic driving force, where they found that the mobility can have a maximum 
at an optimal frequency. More recently, Pizzolato et al. [2] have studied the effects of sinusoidal driving force on 
the polymer barrier crossing over a metastable potential, which is also subject to the reflecting boundary condition. 
They found a similar resonant behavior of the barrier crossing time. In addition, Fiasconaro et al. have studied 
the one-dimensional (ID) polymer chain in the presence of sinusoidal [45] and dichotomically fluctuating [46] driving 
forces. They found that the sinusoidal driving force may induce an oscillating behavior of the translocation time 45j, 
whereas the dichotomic force does not [46] . 

Despite the related work found in the literature, the effect of the polymer-pore interactions on the translocation 
of biopolymers under time-dependent driving forces and in a realistic geometry needs to be studied. In this work, 
the effects of time-dependent driving forces on the translocation dynamics are investigated as a first step towards 
understanding translocation both in vivo and in practical applications. We consider both dichotomically fluctuating 
forces as an example of in vivo noncquilibrium noise |39l 147) and sinusoidal driving forces, which might be easier 
to implement experimentally. We find that the polymer-pore attraction fundamentally changes the behavior of the 
translocation time with respect to the flipping rate of the dichotomic force or the angular frequency of the sinusoidal 
force. For the non-attractive pore, the translocation time has a cross-over to a fast translocation regime at low 
flipping rates (frequencies), but does not have a resonant minimum. For the attractive pore, we show that the system 
exhibits resonant activation within a broad range of physical parameters. We examine the effect of parameters such 
as chain length, driving force and polymer-pore interaction strength on the resonance. The results suggest that in 
vitro experiments with time-dependent driving force might be useful to DNA sequencing. 



II. MODEL AND METHOD 



We consider the translocation of a self-avoiding chain in two dimensions (2D). The polymer chain is mod- 
eled by Lennard- Jones particles interconnected by finitely extensible nonlinear elastic (FENE) springs. Excluded 
volume interaction between monomers is given by the short-range repulsive Lennard- Jones potential: Ut,j(r) — 

4e (f-) 12 — (f) 6 + 6 f° r r < 2 1 / 6 <7 and for r > 2 1 / 6 er. Here, r is the distance between monomers, a is the diam- 
eter of the monomer and e is the depth of the potential well. Neighboring monomers are also connected by FENE 
springs with £/fene(?") = — \kR$ ln(l — r 2 /i?g), where k is the FENE spring constant and R is the maximum allowed 
separation between consecutive monomers. The geometry of the system is shown in Fig. [T] The wall is constructed 
of immobile Lennard- Jones beads of size a. All monomer-wall particle pairs have the same short-range repulsive LJ 
interaction as described above. To investigate the effect of polymer-pore interactions, we consider two main types of 
interactions between the monomers and the pore particles: attractive and non-attractive. In the case of non- at tractive 
interactions, the pore particles are considered to be identical with the wall particles, having a purely repulsive in- 
teraction with the monomers. In the case of the attractive polymer-pore interactions, the cut-off distance of the LJ 
potential between monomer-pore particles is increased to 2.5a (with JTlj constant for r > 2.5a), and the interaction 
strength is characterized by e pm . The interaction can be either attractive or repulsive, depending on the distance of 
the monomer from the pore particles. 

In our simulations, the dynamics of each monomer is described by the Langevin equation 

mr, = -V(C/ L j + C/fene) + F cxt - £v 4 + Ff , (1) 

where m is the monomer mass, £ is the friction coefficient, v, is the monomer velocity and is the random force with 
correlations (F|*(t) ■ Fj*(t')) = 4£kBT5i t jS(t — t'), where ks is the Boltzmann constant and T is the temperature. In 
the pore, the monomers experience an external driving force F oxt = [F + f(t)]x, where F is static (time-independent) 
force, f(t) is the time-dependent force and x is the unit vector along the direction of the pore axis. In this work, 
we consider two types of time-dependent forces f(t). The first is the dichotomic noise, for which f(t) is either 
or —Ad, and changes from one value to the other with flipping rate uj. The dichotomic f(t) has zero mean and is 
exponentially correlated: (f(t)) = and (f(t)f(0)) = A 2 d exp(— 2u)t). As a second example, we consider the sinusoidal 
force given by f(t) = Asm(ujt + (f), where A is the amplitude, lo the angular frequency and <f> is a constant phase. 

We use the LJ parameters e, a and m to fix the scales for energy, length and mass, respectively. The time scale is 
then given by £lj = (ma 2 /e) 1 ^ 2 . The dimensionless parameters in our simulations are Rq = 2, ft — 7, £ = 0.7 and 
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FIG. 1: A schematic representation of the system. The polymer, placed initially on the cis side, is driven through the pore of 
length L = 5 and width W = 3 by the time-dependent external force F + f(t). 



ksT = 1.2. In our model, the bead size corresponds to the Kuhn length of a single-strand DNA, giving approximately 
a sa 1.5 nm. The bead mass is approximately 936 amu, and the interaction strength e corresponds to 3.39 • 10~ 21 J 
at room temperature (295 K). The Lennard- Jones time scale is then 32.1 ps. Where appropriate, we will express our 
results also in terms of the mean translocation time in the absence of f(t), r , which is the physically relevant time 
scale in the system and also easy to measure experimentally. The pore dimensions we set as L = 5 and W — 3, as 
shown in Fig. [I] With the force scale of 2.3 pN, a static driving force of F cxt = 1 then corresponds to a voltage of 
375 mV across the pore (assuming three unit charges per bead and the effective charge 0.094e for a unit charge 48 J. 
The equations of motion are integrated with the Ermak algorithm [351 ISO] with time-step typically At — 0.01, and 
shorter when necessary. 

Initially, the first monomer of the chain is held fixed at the pore (see Fig. [I]) while the remaining monomers are 
allowed to fluctuate until an equilibrium configuration is reached. Then at time t = the first monomer is released 
and the external force is applied. For the dichotomic force, the initial force is randomly selected from +Ad and — 
with equal probability. Correspondingly, for the sinusoidal force, the phase <f> is randomly selected from a uniform 
distribution between [0, 27r]. For small F cxt and weak polymer-pore attraction, the chain may slip out of the pore 
back to the cis side instead of translocating to the trans side. In that case, the equilibration process is repeated and 
the simulation is begun anew. The process is repeated until at least 2000 successful translocation events are recorded. 
In addition to this standard procedure, it is possible to impose a reflecting boundary condition that prevents the first 
bead from slipping back to the cis side. In this case, the simulation is run simply until a successful translocation 
occurs. It turns out that this boundary condition, although widely used in translocation study, fundamentally changes 
the translocation dynamics, as will be discussed in Section [TTT| That is why, unless otherwise indicated, all the results 
presented in this work have been computed without the reflecting boundary condition. 



III. RESULTS AND DISCUSSION 



A. Non-attractive pore, dichotomic driving force 

We begin by considering the purely repulsive polymer-pore interactions, which is the most common case studied in 
the literature. The strength of the Lennard- Jones interaction is e pm = 1 with a cut-off distance of 2 1 / 6 cr. First, we 



consider the dichotomic driving force, with the results for the sinusoidal force presented later in Sec. IIIB We have 
chosen the numerical values F — 0.3 and Ad = 0.2 for the dichotomic force and N = 64 for the chain length, which 
are within the experimental regime. We have checked that within the experimentally relevant force regime and at 
least for N < 128 the qualitative behavior remains the same. 
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FIG. 2: The mean translocation time r and the probabilities Po and P T (see text) as a function of the flipping rate uj of the 
dichotomic force for the repulsive pore. N = 64, F = 0.3, Ad = 0.2 and to ~ 750 ± 4. The statistical error is smaller than the 
symbol size. 



The main results for the dichotomic force as a function of the flipping rate uj are gathered in Fig. [2j As a function 
of the flipping rate uj, we observe two distinct regions. In the fast flipping regime, uj 3> 1/to, the average translocation 
time is t{uj) ~ tq. Here, To is the translocation time in the absence of dichotomic forces, i.e., Ad = 0. In this limit, 
due to the high flipping rate, f(t) changes its sign many times during the course of the translocation and is averaged 
out to zero over the whole process. Therefore, we have t(uj) rs tq for uj — > oo. In principle, in this limit, the time- 
dependent force becomes a rapidly fluctuating (5-correlated noise similar to the thermal random force Ff-. For the 
monomers inside the pore, the modified correlation of the random force is given as ([f(t)x + F R (t)} ■ [/(0)i + F R (0)]} = 
(4£fc^T + A 2 d /uj)5(t) = A^ksTdif), where last expression is obtained in the limit uj — ► oo. Therefore, the effect of the 
dichotomic force in this limit is vanishing, and we recover t(uj) w tq, as shown in Fig. [2] This result is also in agreement 
with Refs. [4"4ll46| . As the flipping becomes slower, uj < 1/tq, we observe a cross-over to a faster translocation regime, 
with t{uj) < tq. This result is in sharp contrast with Refs. |44j446j . where it was found that r(w) > To. In addition, 
we do not find a global minimum of t(uj) at any finite u, unlike Refs. [4"4l I45j . 

To understand the behavior of t(uj) at small u, we need to look at the probability of achieving a successful 
translocation. Due to confinement within the pore, the chain experiences an entropic free energy barrier |12l I13L 115) , 
as illustrated in Fig. [3] Because of fluctuations, there is a finite probability that the chain slips back to the cis side 
instead of translocating to the trans side. Therefore, the probability of translocation is less than one and increases 
with increasing driving force (for details of the translocation probability as a function of various system parameters, 
see Ref. [22]). Thus, within the set of successful translocations, we expect to find a larger number events that have 
positive fit), as compared to those with negative f(t). We characterize this dependence of the translocation probability 
on the driving force by looking at the set of successful translocations, from which we calculate the distribution of 
f(t) at the beginning of translocation (t = 0) and at the final moment of translocation (t = t). The probabilities 
Po = P[f(0) > 0] and P T = P[/(r) > 0] that the force f(t) is positive for t = and t = t, respectively, are shown 
in Fig. [2j In the high-rate regime, the flipping rate is too high for f(0) or /(t) to be correlated with the chain 
dynamics, and therefore Po and P T approach 0.5. On the other hand, in the low-rate regime, the positive direction 
of /(0) is strongly favored (Po ~ 0.98). In addition, since the correlation time of the driving force is much longer 
than r, the driving force remains constant during the whole translocation process with high probability, being either 
-fcxt = F + Ad or P cxt = F — Ad- Therefore, in this limit, the average translocation time is given by the weighted 
average 

t = P t + + (1-P )t_. (2) 

Here, r + and r_ are the translocation times with the total force F + Ad and F — Ad, respectively. Assuming that the 
translocation time is inversely proportional to the driving force, r(/) ~ Eq. ^ gives t(uj) s» 0.65t in the low-w 
limit. This agrees well with the results in Fig. [2j Therefore, the cross-over to the fast translocation regime (r < To) 
as the flipping becomes slower is simply explained by the fact that for low flipping rate the chain is most likely to 
translocate when /(0) > 0. This strong bias for selecting the initial value /(0) induced by the entropic barrier is the 
crucial difference between this work and Refs. [H 46 In Refs. |4"4"144l)] this kind of selection does not occur because 
the translocation probability is one, independent of the time-dependent driving force. A similar effect can be obtained 
in our model by imposing a reflecting boundary condition that prevents the first monomer from exiting the pore to 
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FIG. 3: Illustration of the free energy of the polymer chain as a function of the number of translocated monomers s. The 
dotted line indicates the free energy for the non-attractive pore, which has no well structure. A reflecting boundary condition 
at s — forms a free-energy well (blue shaded area). Attractive polymer-pore interactions can also create a free-energy well 
(schematically shaded red). 




FIG. 4: The distribution of translocation times for chain length N — 64 and F = 0.3 under dichotomic driving force in the 
non-attractive pore. Panel a) shows the distribution for Ad = 0, while panels b) - d) show the distribution for Ad = 0.2. 

the cis side. However, we stress that this kind of boundary condition may not be realistic for, e.g., the translocation 
of a ss-DNA molecule through a pore, although it has been used in many studies. 

Finally, we look at the distribution of translocation times. In the high flipping rate limit, the distribution is very 
similar to the zero amplitude case (see. Fig. [4]). In general, the distribution at this limit is either Gaussian (for 
large enough F) or has an exponentially decaying tail (for small F). In the present case, the distribution is almost 
Gaussian with a slightly elongated tail, which differs greatly from the typical distributions of thermally activated 
processes. Furthermore, as shown in Fig. [3j there is no metastable well (pretransition state) within which the chain 
attempting to escape would oscillate. Thus, in this case, the resonant minimum of t(uj) does not exist, in contrast 
to Ref. [44], where the adopted external potential has a pretransitional well. At lower flipping rates, the peak of the 
distribution moves toward shorter translocation times, as the trajectories having f(t) predominantly in the positive 
direction are favored (signaled by increasing Pq and P T ). However, also the probability of long translocation times 
increases. These events correspond to the trajectories with negative f(t). As the flipping rate is further decreased, 
most of the successful translocations occur with f(t) > 0. In the low-rate limit, one retains two peaks, corresponding 
to F ext = F + Act and F cxt = F — For F = 0.3 and A^ — 0.2, only the former is practically visible. 

B. Non-attractive pore, periodic driving force 

As a second case, we study the translocation through a non-attractive pore under sinusoidally time- dependent 
driving force f(t) — Asm(u)t + 4>). The average translocation time t(uj) and the probabilities Pq and P T are shown in 
Fig. [5] as a function of the angular frequency ui. For comparison with the dichotomic case, we use the parameter values 
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FIG. 5: The mean translocation time r and the probabilities Po and P T as a function of the angular frequency uj for the periodic 
force and repulsive pore. N = 64, F = A = 0.3 and to ~ 750 ± 4. The statistical error is smaller than the symbol size. 



N = 64, F = 0.3 and A — 0.3. The time-averaged amplitude of the time-dependent force is then (|j4sin(wi)|) t = 
2A/tt s» 0.2, which corresponds to the value of Ad used in the previous Section. In the low-frequency (w <§; 1/to) and 
high-frequency (w 3> l/ T o) limits we obtain results similar to the dichotomic force explained above: in the high-w 
limit, t(uj) ro and, in the opposite limit of small u>, t(uj) < tq. The average translocation time is given by a relation 
analogous to Eq. 

r = i / p(0)r(0)d0. (3) 
27T Jo 

Here p((f>) is the probability density of the initial phase cf> within the set of successful translocations and r((/>) is the 
average translocation time corresponding to the driving force F + Asin(<p). Similarly to the dichotomic case, the 
distribution p(<j)) is uniform only in the high frequency limit, while in the low frequency limit, values of <j) giving 
/(0) > are strongly favored (cf. Fig. [5]), as we shall see. This leads to larger average driving forces and consequently 
faster translocation. 

In the intermediate regime (ui ~ 1/tq), the periodic time-dependence of the driving force fundamentally affects the 
translocation dynamics. Instead of a simple cross-over in t(u), one gets a series of local minima and maxima. In 
addition, the probability P T has a local maximum in the vicinity of a local minimum of r. In many cases, these could 
be argued to indicate resonant activation 40 . However, in this case they have a deterministic origin weighted with the 
distribution p(4>). To show this, we consider a coarse-grained model for the translocated segments already studied in 
Ref. |12[|13j . with the entropic contributions therein neglected to make it analytically solvable. This approximation is 
reasonable because apart from the short initial (and final) stage of translocation, the entropic force is small compared 
to the mean driving force F. Our model is the ID equation of motion for (s(t)), the average number of translocated 
segments, under the sinusoidal driving force with a fixed value of </>, 

Ceff^^=F[l+8in( W t + ^] ) (4) 

where £ e ff is the effective friction. Eq. Q can be analytically solved for (s(t)) with the initial condition s(0) = 0. 
Specifically, we are interested in the time that it takes for the system to evolve from s = to s = N as a function 
of the phase, t(4>). We fix £ e g by setting the time-scale of the model so that F/£ c s — N/tq, giving £ e ff w 17.6 for 
N = 64 and F = A = 0.3. The integration of Eq. Q yields 

T (4>) = T o + — [cos(u;t + (f>) - cos(0)] . (5) 

This describes the approach of the translocation time r to ro in the u — >■ oo limit, as well as the local oscillation in 
the intermediate u> regimes. Once t(</>) is obtained as a function of <j> as well as u from Eq. |5]), the translocation time 
averaged over <fr is found from Eq. © 

In Fig. U we compare our model with the JV-particle Langevin dynamics simulations. First, the dotted line shows 
the results for a uniformly distributed <f>. In contrast to the LD simulations, this curve shows a global minimum of 
translocation time, and also a strong oscillating behavior as a function of uj. The behavior is very similar to the simple 
ID chain model driven by sinusoidal force studied in Ref. 45J. 
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FIG. 6: Comparison between LD simulations (N = 64, F — A — 0.3) and the theoretical toy model. Dotted line: toy model 
with uniformly distributed <f>, solid line: toy model with Boltzmann distributed <j> (see text). The latter shows good agreement 
with the LD results (circles). 

However, <fi should not be chosen uniformly. In the properly formulated translocation problem, the chain has to 
overcome the initial free energy barrier, which leads to nonuniform distribution of 4>. In the zero-frequency limit, the 
translocation probability follows the Boltzmann distribution, which depends exponentially on the height of the initial 
free energy barrier (cf. Fig. [3]). Hence, we put the distribution in the form p(4>) ~ exp [asin(c/>)]. Since we consider 
only the processes that complete the translocation, a is a nontrivial function of not only fc^T, F, A but also u). In 
our procedure, the parameter a is obtained by fitting the integral J* p(4>)d<j> to the probability Po(uj) for each uj. The 
a, obtained as a(oj) = (15u> + 1/2.6) -1 , serves as an empirical interpolation between the Boltzmann distribution for 
u) <C 1/to and the uniform distribution of 4> for iu 3> 1/to- This a then gives the distribution p(4>), over which the 
average of t(c/>) is taken by Monte Carlo integration to eventually find the average translocation time r. The results of 
our model with this distribution of are shown in Fig. [6] as a solid curve. The model reproduces the essential features 
of the full A^-particle Langevin dynamics simulation: the cross-over to fast translocation as uj is decreased, and the 
global and local maxima of r near ui ~ k/tq. This exercise clearly shows that the local maxima and minima are a 
result of deterministic dynamics and the nonuniform distribution of c/>, and are not indications of resonant activation. 

The difference between the sinusoidal and dichotomic driving forces can also be identified in the translocation time 
distributions. In the high and low frequency limits, one recovers distributions very similar to the dichotomic case. In 
the intermediate frequency regime, on the other hand, the sinusoidal time-dependence shows as a periodic modulation 
of the underlying distribution. Here, the distribution has multiple peaks, which correspond to translocations occurring 
when f(t) > with high probability. Each peak corresponds to one period Tq = 2tt/uj of the sinusoidal force, with the 
distance between consecutive peaks being At/tq w 2tt/lut$, as shown in Fig. [7](b). Near ui ~ vr/ro, where the average 
translocation time achieves its maximum, the distribution shows two distinct peaks, corresponding to fast (r < To) 
and slow (r > t ) translocation (see. Fig[7]d). The leftmost peak corresponds to events that occur roughly between 
Tq/4 < t < To/2, with a phase <j> between < <j> < ir/2. For these events, f(t) is positive for the whole translocation 
process, resulting in faster than average translocation. The peak on the right, on the other hand, corresponds to the 
events with 3Xh/4 < r <Tq and 7r/2 < < n. Here, although f(t) starts positive, it quickly crosses over to negative 
values. Typically, translocation occurs when f(t) turns back to positive. Thus, the average f(t) during one event is 
negative, giving longer than average translocation time. As the frequency uj is decreased, the rightmost peak becomes 
smaller as the phases <f) corresponding to that peak become less probable. As a result, the average translocation time 
crosses over to the regime where r < tq. 

C. Attractive pore, dichotomic driving force 

We have shown above that for purely repulsive pore-monomer interactions, the system does not exhibit resonant 
activation. This is due to the absence of a proper free-energy well, in which an attempt frequency of crossing the 
imminent barrier is well defined. Introducing attractive interactions between the polymer and the pore modifies the 
free energy in such a way that a well is formed (schematically shown in Fig. [3| , and translocation becomes a thermally 
activated barrier crossing process [TH1 H21 US] . Therefore, for the attractive pore, we expect to find a resonance similar 
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FIG. 7: The distribution of translocation times for chain length N — 64 and F — 0.3 under sinusoidal driving force in the 
non-attractive pore. Panel a) shows the distribution for A = 0, while panels b) - f) show the distribution for A = 0.3. 
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FIG. 8: The mean translocation time r and the probabilities Po and P T for the dichotomic force and attractive pore. N = 32, 
F — 0.5, Ad = 0.2, e pm = 1, and tq « 226.8 ± 0.6. The statistical error is smaller than the symbol size. 

to that reported for the polymer escape in Ref. We start with the case of dichotomic driving force, which is 

somewhat more pedagogical than the sinusoidal force case. 



1. Dependence on the polymer-pore interaction strength e F 



First, we study the effect of the polymer-pore interaction strength e pm on the average translocation time r. For the 
attractive pore, we use the value 2.5a for the cut-off distance of the Lennard- Jones potential, which yields an attractive 
force between the pore and the monomer at distances 2 1 / 6 er < r < 2.5a. The chain is driven by a dichotomically 
fluctuating force with the flipping rate uj and correlations described in Section [TTJ In Figs. [8] and |9j we show the 
average translocation times t(uj) for the chain length N = 32, with the polymer-pore interaction strength e pm = 1 
and e pm = 2.5, respectively. Here, the static force F = 0.5 and the amplitude of dichotomic force is Ad = 0.2. In the 
high flipping rate regime, U) 3> 1/tq, t(cj) ~ To, as for the non-attractive pore. On the other hand, for u) -C 1/tq, the 
translocation time is r > To. This behavior is completely opposite to the non-attractive pore case. Nevertheless, it 
can be explained by the same arguments. The average translocation time is given by Eq. Q. However, for sufficiently 
large e pm the selectivity with respect to the initial driving force /(0) is fairly weak, because a strong attraction between 
the pore and the polymer prevents the escape to the cis side. For example, for e pm = 1.0, Po ~ 0.63, as shown in 
Fig. [8] Assuming inverse dependence of the translocation time on the driving force, Eq. ([2| gives r 1.06to, which 
is in agreement with the simulation results. 
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FIG. 9: The mean translocation time r and the probabilities Po and P T for the dichotomic force and attractive pore. N = 32, 
F — 0.5, Ad = 0.2, e pm = 2.5, and to ~ 2202 ± 29. The statistical error is smaller than the symbol size. 



In the intermediate regime (w rs \/tq), the translocation time t(oj) shows different behavior depending on the 



value of e 



pm ■ 



While for e 



pm 



1. 



t(u>) monotonically decreases as ui increases, for e pm = 2.5, t(w) has a minimum 
at an optimal flipping rate ujtq w 1.8. Related to this, we obtain the probabilities Pq(ui) and P T (w). For e pm = 1, 
Po monotonically increases as ui decreases, similarly to the non-attractive case, but only by approximately 0.1. For 
e pm = 2.5, Po f=a 0.52, almost independent of ui. On the other hand, P T shows nonmonotonic behavior, having a 
maximum at ujt w 0.4 and ujtq rs 1.0 for e pm = 1 and 2.5, respectively. Typically, for a barrier crossing problem, such 
a maximum is an indication of resonant activation, and is accompanied by a minimum in the crossing time |401 147) . 
However, out of the two cases, e pm = 1 and e pm = 2.5, only in the latter has a minimum in t(cj). Furthermore, for 



e pm = 2.5, the flipping rates uj at the maximum of P T {u) and at the minimum of t(uj) do not coincide. To understand 
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these results, we divide the translocation process into three components [HI [Ml [2S] : ^ initial filling of the pore, 2) 
transfer of the polymer from cis to trans side, and 3) emptying of the pore, as shown in Fig, 
time is then r = t± + t 2 + r 3 , where 73 is the time for the ith process. In Fig. 11 
e pm = 1 and e pm = 2.5. For the larger e pm , t 3 dominates the translocation time 



The translocation 
r i + T2) and T3 for 



we show Ti,2v 
. As ui increases, t 12 {ui) decreases 



gradually, but ra(cj) behaves non- monotonically. In addition, the minimum of T3 coincides with the maximum of 
P T (for e pm = 1 this is barely observable). This indicates that P T {u) and T3(u>) are highly correlated. Thus, the 
nonmonotonic behavior of t(ui) occurs because the time-dependent force couples to the pore emptying process, i.e., 
the crossing of the final free-energy barrier (cf. Fig. [3]). The coupling to the first two processes is very weak, and 
does not significantly contribute to the resonant activation. However, since slightly decreases as ui increases, 
the optimal flipping rate that yields the minimum of translocation time r is somewhat larger than the rate at the 
minimum of T3 . 

The translocation time distribution P(r) also profoundly depends on the magnitude of e pm , as shown in Fig. 12 The 
left column shows the case with e pm = 1, while the right column corresponds to e pm = 2.5. The first row shows P(t) 
in the presence of static force F only (corresponding to very high flipping rate uj). While for e pm = 1, P(t) is nearly 
Gaussian centered at To, for e pm = 2.5 the distribution is an exponential. This indicates that strong polymer-pore 
interactions make translocation an activated process, where the chain has to surmount the final free-energy barrier 
before it can completely translocate to the trans side (cf. Fig. [3]). For e pm = 1, as ui decreases, P(r) gradually splits 



into two Gaussian distributions, centered at 



-tq and F _ . t . For e pm = 2.5, as uj decreases, P(r) is changed in a 



F+A d ' u F-A"->- c P m 

nontrivial way: at intermediate flipping rate uj w 1/to, the tail of P(r) is shortened, but for lower ui, P(r) develops a 
long tail. The behavior of P(t) at intermediate uj is closely related to the probability P T {u>) in Fig. [9} Although /(0) 
is either positive or negative with similar probability, most of the successful translocations finish with /(r) = +Ad, 
which results in a shorter translocation time. This is the reason for the small probability of long translocation times. 
On the other hand, at very low ui, P(t) becomes a combination of two exponential distributions, each corresponding 
to the translocation time with the driving force either F + Ad or F — Ad, which results in sharp increase of t(uj). 



2. Dependence on the chain length N 



In the previous section, we found that for large e pm , the translocation time r has a resonant minimum as a 
function of the flipping rate ui. Here we study how this behavior changes with chain length N. Figure [13] shows 
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FIG. 10: The translocation process divided into three stages: 1) initial filling of the pore, 2) transfer of the polymer from the 
cis side to the trans side, 3) the final emptying of the pore. The corresponding times of the subprocesses are n, T2 and T3, with 
the total translocation time r = T\ + r 2 + r 3 . 




FIG. 11: The times 7-1,2 (dashed lines) and T3 (solid lines) for e pm = 1 (squares) and e pm = 2.5 (circles). Here N = 32, F — 0.5, 
Ad — 0.2 for both cases. While n,2 monotonically decreases as uj increases, T3 shows a resonant minimum. The statistical error 
is smaller than the symbol size. 



the translocation times t(u>) for 16 < N < 128, with e pm = 2.5, F = 0.5, and Ad = 0.2. The optimal flipping 



rate that yields the minimum of translocation time is roughly independent of N (see the inset of Fig. 13). Since 
the resonant behavior occurs during the last emptying process, this indicates that T3 is independent or only weakly 
depends on the chain length. The free energy barrier of the last emptying process can be approximated as A_F — 
L(e pm — F/2 — g(N)) [HI [22J [25] . Here the first term accounts for the polymer-pore interactions, the second term for 
the potential energy difference across the membrane due to the driving force, and the last term is due to the entropic 
free energy [TU [T31 [T5]. For the pore-emptying process the entropic force g(N) is in the positive direction, slowly 
increasing with N and eventually saturating for very long chains On the other hand, as shown in Ref. for 
the static driving force and N < 200, ti 2 approximately increases as ~ N 1 - 5 . The normalized translocation time 
is 

t{uj) _ T2M + t 3 (w) _ t 1i2 (u)/t 3 (0) + t 3 (uj)/t 3 (0) 
t ~ Ti l2 (0)+7a(0) n )2 (0)/r 3 (0) + l ' U 
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FIG. 12: The distribution of translocation times for the dichotomic force and attractive pore. N = 32, F = 0.5, and Ad = 0.2. 
The left column shows the distributions for e pm = 1 and the right column for e pm = 2.5. 
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FIG. 13: Translocation times for chain lengths 16 < N < 128 with the dichotomic force and attractive pore. F — 0.5, Ad = 0.2 
and e pm = 2.5. While the optimal rescaled flipping rate (ojto) shows a slight dependence on N (main figure), the unnormalized 
flipping rate (w) is independent of N (inset). The statistical error is smaller than the symbol size. 



where Ti(0) is time for the ith process in the absence of the time-dependent driving force f(t). In the short chain 
limit, Ti 2,(w) <C t 3 (oj) : so that the normalized translocation time is ' r ^- ~ ^(o) ' ^ n ^ e °t ner hand, in the long 
chain limit, ti*s(u) ^> T3(u>), giving ps ■ These limiting situations predict that for short chains, one should 

observe a strong minimum in t(cj), whereas for very long chains, the minimum should vanish. This trend can be 
observed in Fig. |13[ where the minimum of translocation time becomes less pronounced as iV increases. In addition, 
the optimal flipping rates are quite independent of N . This is in contrast to the results of Ref. jSHj; where the authors 
consider the translocation of a rigid rod in the presence of a reflecting boundary condition at s = (see Fig. pjl). 
In that case, all the segments of the polymer are subject to the external forces, making the translocation time very 
sensitive to their minute changes and the optimal flipping rate decreases with N . However, in the present case the 
number of segments within the attractive pore remains small throughout the translocation process. Thus, the effect 
of the external forces becomes small as the chain gets longer. 
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FIG. 14: Translocation times for driving forces 0.5 < F < 4 for the dichotomic force and attractive pore. Ad = OAF, e pm = 2.5. 
In this case, the optimal rescaled flipping rate (wto) is roughly independent of F (main figure), while the unnormalized flipping 
rate (lu) strongly increases with F (inset). The statistical error is smaller than the symbol size. 



3. Dependence on the driving force F 

Next, we study the effect of changing the driving force magnitude. We consider static driving forces between 
0.5 < F < 4 with the amplitude Ad fixed as Ad = OAF. The results are shown in Fig. [Mj One can see that the 
optimal flipping rate for the resonant minimum translocation time increases with increasing F. In addition, the 
resonant minimum becomes shallower, finally disappearing for large F and t(w) becomes a monotonic function of u>. 
It is because the free energy barrier of the last emptying process vanishes for large F. It is of interest to study the 
critical driving force F c , for which r(w) changes from non-monotonic to monotonic. F c can be approximated from 
the condition AF — 0, so that F c — 2(e pm — g(N)). For N — 32 and F = 0.5, t{uj) becomes monotonic for e pm ^ 1 
(cf. Fig. [8]), giving the estimate g(N) s» 1. Therefore, F c w 3 for e pm = 2.5 and N — 32. This estimate seems to 
be reasonable as shown in Fig. [14] This result also shows that the condition of the non-monotonic behavior of the 
translocation time is determined by the competition of the polymer-pore interaction e pm and the driving force F. 



4- Dependence on the driving force amplitude Ad 

As the last case of the dichotomic force, we study the effect of changing the dichotomic force amplitude Ad while 
keeping the static driving force F constant. The results for the translocation time are shown in Fig. [15] With 
increasing Ad, the resonant minimum becomes deeper and the resonance flipping rate uj gradually increases. For very 
large Ad, the resonance disappears and the translocation time becomes a monotonic function of w, similar to the 
non-attractive pore case. This shift in behavior is because the selectivity of initial sign of /(0) becomes stronger for 
larger Ad- The transition to this regime happens when the initial barrier (see Fig. [3]) that prevents the chain escape to 
the cis side becomes comparable to the thermal energy and the escapes become frequent. For the negative dichotomic 
force, /(0) = — Ad, the barrier can be written in the form AF C j S = e pm + (F — Ad)/ 2 — g(N). The pore length L does 
not enter the relation because in the initial configuration, only the first bead is inside the pore. For e pm = 2.5 and 
ksT = 1.2, the requirement Ai 7 ^ rs ksT gives the estimate Ad ~ 1 for the transition from the nonmonotonic t(cj) 
to the monotonic one. This estimate matches the data in Fig. [T5| 



D. Attractive pore, periodic driving force 



For the attractive pore, the dichotomic and periodic driving forces give qualitatively very similar results. Also 
for the periodic driving force, f(t) — Asin(ujt + </>), a resonant minimum of the average translocation time appears, 
in contrast to the non-attractive pore case. In addition, the location (frequency u) and the depth of the minimum 
depends on the various parameters in essentially the same way as described in the previous section, which indicates 
that the origin of the resonance is the same: the time-dependent force being most co-operative to translocation during 
the pore-emptying time T3. However, there are also some obvious differences. For the periodic driving force, oscillatory 
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FIG. 15: Translocation times for the dichotomic force and attractive pore for amplitudes Ad G {0.2,0.4,0.8,1.6}. F = 0.5, 
e pm = 2.5. The statistical error is smaller than the symbol size. 




FIG. 16: Translocation time t as a function of frequency ui for the periodic driving force f(t) = A sin(ujt+<j)) for 0.1 < e pm < 2.5. 
F — 0.5, A — 0.3 and N = 32. The inset shows a magnification of the data for 0.1 < e pm < 1.5. The statistical error is smaller 
than the symbol size. 



behavior similar to the one described in Section IIIB emerges, in addition to the resonant activation. In this Section, 
we will briefly describe the essential differences between the two driving schemes and examine some of the implications 
of employing the periodic driving force. 



1. Dependence on the polymer-pore interaction strength e F 



First, it is instructive to consider the dependence of the translocation time t(u>) on the strength of the polymer-pore 
interaction strength e pm . As shown in Fig. 16 for low interaction strengths, one recovers the transition from fast to 



slow translocation with local minima and maxima in t(u>), characteristic of the non-attractive pore case. For larger 
e pm , a resonant minimum in t{uS) develops, similarly to the dichotomic driving force. However, the global maximum 
observed for the non-attractive pore persists, although it is reduced to a local maximum located within the resonance 
minimum. This local maximum arises because of the interplay of the periodic forcing and the non-uniform distribution 
of the phase 4>, as we will discuss below. 
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FIG. 17: Translocation time r as a function of frequency u for the periodic driving force f(t) — Asin(ujt + (f>) for amplitudes 
A <E {0.3,0.6,0.8,1.0,1.2,1.5,1.8,2.4,3.0}. Other parameters are F = 0.5, e pm = 2.0 and N = 32. The inset shows the 
dependence of the frequency ai m in of the global minimum translocation time on the amplitude A. Here r w 500 ± 6. The 
statistical error is smaller than the symbol size. 



2. Dependence on the driving force amplitude A 



To highlight the differences between the sinusoidal and dichotomic driving forces, we look at how the translocation 
time r(w) changes with the driving force amplitude A. For the dichotomic force, as Ad is increased, one merely crosses 
from the non-monotonic t(w) with the resonant minimum to the monotonic t(lS) characteristic of the non- at tractive 
pore case. The sinusoidal driving force, on the other hand, exhibits much richer behavior. As shown in Fig. |17| as 
the amplitude A is increased, the minimum becomes deeper and slowly moves toward higher frequencies. In addition 
to the original one, another (local) minimum appears at the low-frequency end of the spectrum and travels down the 
t(uS) curve as A is increased. Eventually, the new minimum becomes a global one. This produces a sudden transition 
in the frequency of minimum translocation time, w m i n , as shown in the inset of Fig. |17| Finally, at sufficiently large 
A, the new minimum merges with the original one. To better understand this complex behavior, we again divide the 
translocation time r to the three components t±, T2 and T3. Looking at T\ t 2 = T\ + T2 and T3 separately reveals that 
the original global minimum of t(uj) is associated with the pore emptying time T3, as shown in Fig. 18 The additional 
minimum, on the other hand, is related to the periodic back-and-forth movement of the chain, which is visible as a 
nonmonotonic behavior in t\^{ui). 

Let us first examine the pore emptying time T3, since that shares many similarities with the dichotomic force case. 
In both cases, the minimum of t^(lo) occurs due to resonant activation. At the corresponding resonant frequency, 



the probability P T also reaches a maximum (not shown), similarly to the dichotomic force case. As shown in Fig. 18 



for the sinusoidal force, there is also a small local maximum in T3(oj) at ujtq « 15. Surprisingly, here P T also has a 
local maximum, which should indicate efficient crossing of the final barrier. However, instead of the expected decrease 
in T3, one sees a slight increase. The reason is that there is a special mismatch between the frequency oj and the 
translocation time so that the period Tq = 2it/uj s» T1.2. In other words, typically it takes the chain one period of 
f(t) just to traverse from its initial position to the configuration where it may try to surmount the final free-energy 
barrier (cf. Fig. [3]). Thus, the chain essentially misses the first opportunity to cross the final barrier, which slightly 
increases T3. This can be also seen as a suppressed first peak in the translocation time distribution of Fig. 19 b) . 

Finally, let us look at the translocation times T\ and T2. The combined time t\^% shows features similar to the 
non-attractive pore case, where the periodic driving force produces a series of alternating minima and maxima. In 
the case of the attractive pore, the selection over the initial phase is weaker because of the free-energy barrier that 
prevents the escape to the cis side. Consequently, the qualitative behavior of the ti^{oj) curve is closer to the model 
of Eq. Q with uniformly distributed phase (see Fig. j6j). Essentially, the local maximum in n,2 is produced by the 
interplay of w-dependence of the distribution of cf) and the periodicity of the driving force. Close to the resonant 
minimum, u>t q m 5, the distribution is bimodal, as shown in Fig. |19| Similarly to the non-attractive pore case, the 
first peak corresponds to the events that occur within the first half-period of the force fit), i.e., < r < Tq/2, 
and whose initial phase is typically — 7r/4 < <j> < ir/2. Therefore, for these trajectories, f(t) > for most of the 
process, and translocation occurs faster than average. In contrast, the second peak corresponds to tt/2 < <f> < 7tt/A 
and Tq/2 < r < To. As the frequency u decreases, the second peak moves further to the right (Fig. [l|fcl),(e)). 
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FIG. 18: The translocation time r and its components n,2 and T3 as a function of uj, showing that the leftmost minimum in 
t(uj) is associated with ti,2- Parameter values used are F — 0.5, A = 1.0, e pm = 2.0 and N = 32. The statistical error is 
smaller than the symbol size. 
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FIG. 19: Distribution of translocation times for the periodic driving force with N = 32, e pm = 2.0, F — 0.5, A — 0.0 (panel a) 
A = 1.0 (panels b-f). 



This increases the average translocation time. At the same time, the area under the peak decreases, because with 
decreasing uj, the distribution of <fi becomes less uniform, favoring <fi belonging to the first peak. This tends to decrease 
t. The combination of these two factors creates the maximum of n^. For larger A, the selection over 4> is stronger, 
so the second factor starts to dominate already at relatively high frequencies. Conversely, for small A, the second 
peak in P(t) persists for even very small uj. This explains why the maximum occurs at lower frequencies for small A, 
and moves towards the high-frequency end as A is increased. 



IV. SUMMARY 



In this work, we have studied the translocation of polymers under a time-dependent driving force using Langevin 
dynamics simulations. In particular, we have extracted the dependence of the average translocation time on the 
flipping rate ui of the dichotomic driving force and the corresponding dependence on the angular frequency ui for the 
sinusoidal driving force. We have also examined the influence of various other physical parameters on the translocation 
dynamics. 

We find that the interactions between the polymer and the pore play a fundamental role in the dynamics of the 
system. For the non-attractive interactions, the translocation time shows a cross-over to a faster translocation regime 
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at low flipping rates. For the sinusoidal force, in addition to the cross-over, we observe a series of local minima 
and maxima, produced by the periodicity of the driving signal. However, in this case we do not observe a global 
minimum of the translocation time for any finite uj. On the other hand, with attractive polymer-pore interactions, 
which represent naturally occurring biological pores such as the a-hemolysin, the situation is very different. In this 
case, the translocation becomes a thermally activated process due to the attraction between the pore and the polymer. 
Optimal modulation of time-dependence driving force induces a resonant activation, manifesting as a global minimum 
in the translocation time at finite uj. We also find that, although the details of this resonance depend on various 
system parameters, in general the resonance is quite robust and occurs for both the dichotomic and sinusoidal driving 
force. Typically the resonant flipping rate (angular frequency) is found in the neighborhood of w ~ 1/to, with To being 
the translocation time without the time-dependent component of the driving force. For an experimentally typical 
translocation time of the order of 100 /is [5], this corresponds to the rate (frequency) in the kilohertz regime. 

Theoretically, the occurrence of the resonance relies on the existence of a free energy well, from which the poly- 
mer escapes via thermal activation. In practice, to observe the resonant behavior, one should choose the physical 
parameters so that one has: 1) strong enough polymer-pore interactions, 2) relatively short chain length, 3) small 
enough static driving force. For example, a poly(dA) 10 o chain driven through an a-hemolysin pore should display the 
resonance for pore voltages of roughly V 1000 mV. For V ~ 200 mV, we would expect to find the resonance in the 
neighborhood of uj ~ 1 — 10 kHz at room temperature. Furthermore, in the case of the dichotomic driving force, one 
should also have a relatively small amplitude of the time-dependent force, whereas for the sinusoidal force, even sig- 
nificantly larger amplitudes can still produce the resonance. In the latter case, a more complicated behavior emerges, 
as the driving force not only assists translocation during the pore-emptying time T3, but also significantly alters to 
the overall motion of the chain. Our findings suggest that time-dependent driving forces may play a fundamental part 
in polymer translocation in biological systems, and may also be useful in practical applications such as sorting and 
sequencing of DNA molecules. 
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